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Abstract 

(N , 

^ We have made the comparison of three parametrizations of the proton electromagnetic 

form factors in the space-Uke region using the largest data set on e^p elastic scattering 
ever used. All models have a correct analytical structure with a minimal set of well known 
singularities. For the first time an effective two-pion intermediate coulombic bound state 
I (pionium) was included into GVMD for nucleon formfactors. For the proton charge radius 

(N 



X 



an estimate derived from our best fits is as follows: 

r|[fm] = 0.897 ± 0.002(exp) ± O.OOl(norm) ± 0.003(models). 



I The traditional "model independent" fit by quadratic polynomial gives, in our estimation 

■ technology, a much cruder estimate 



r 



p 

E .polinomial 



[fm] = 0.887 ± 0.006(exp) ± 0.096(norm). 



Oh! 

(-H \ Both estimates are closer to the value estimated from the recent hydrogen Lamb shift mea- 

surements. 



Introduction 

Since the discovery of the deviation of differential cross sections for ep elastic scattering |]l|] from 
that of calculated for point-like nucleons, numerous efforts were devoted to measurement of 
the "spatial extention" of nucleons and the corresponding electromagnetic form factors (FFs). 
There are more than 50 experimental reports with data on e^p elastic scattering, starting from 
a pioneering study of Hofstadter and co-workers in 1953-1956, up to recently reported new 
experimental data of Andivahis et al. in 1994 |]lO| (see Table 5.1). 

Many attempts were made to model charge and magnetic moment distributions in the nu- 
cleons, numerous phenomenological parametrizations of (FFs) were proposed. The majority of 
proposed parametrizations is based on various modifications of generalized vector mesons dom- 
inance (GVMD) model guided by the requirements from general principles of quantum field 
theories and the modest amount is motivated by QCD features. The common features of the 
parameterizations are the requirements of providing reasonable values of the rms "charge" and 
"magnetic" radii of the nucleon (t ^ asymptotic ), as well as reproducing QCD-like behaviour 
at large \t\. 
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Unfortunately, it is still impossible to derive the FFs behaviour, even in the leading asymp- 
totic terms, from the "first principles" and any extraction of derived physical parameters (such 
as radii and constants in an asymptotic behaviour at large \t\) would be "model dependent." 

On the other hand, the precise knowledge of the nucleon FFs with as accurate as possible 
errors treatment is needed in a number of phenomenological applications and even in precise 
tests of well established theories. The analyses of charged hadrons elastic scattering in the 
Coulomb-nuclear interference regions and of the data on Lamb shifts in atomic hydrogen are the 
well known examples of such applications. 

It is often argued that the extraction of FFs from the experimental data is an ill posed 
problem. The phenomenological basis of experimental FFs extraction is the well known Rosen- 



bluth formula for da/dt in the one photon exchange approximation |11]. Different experimental 
groups have extracted FFs from their own data on elastic scattering with loose information on 
the sources and on the nature of different errors estimated for differential cross sections. As a 
result we have no reference numbers to make all the data on FFs extracted by different groups 
mutually consistent. 

As the common practice of usage the experimental information on FFs in applications is via 
the "currently best" parametrizations we propose in this paper a schema to select the best FFs 
parametrizations by fitting to the world e^p elastic scattering data without using the derived 
data on FFs. It makes possible to have a relatively complete control of the error matrices of 
the technical and physical parameters and to estimate the errors induced by switching of the 
competitive parametrizations and errors induced by our method used to make world database 
self consistent in relative normalizations of different experiments. 

In this paper we analyze only e^p elastic scattering data. The article is organized as follows. 

Section 1. Here we reproduce a compilation of traditional formulas used in data analyses 
for extraction of the FFs from the data on da/dt. 

Section 2. presents the description of our preparatory handling of the experimental data 



published so far and stored in the ReacData database of PPDS system. Our procedure of 



the cross assessments of data and models is also outlined in this section and in the Appendices. 

Section 3. is a compilation of the parameterizations used in cross assessments. 

Section 4. Here we describe the preparatory fit results for each model selected for cross 
assessments. 

Section 5. gives details of our models driven procedures to make data from different exper- 
iments consistent with each other. 

Section 6. Here we describe the procedure of the model-dependent filtration and present 
the final results of the fits to the filtered data. 

Section 7. The summary of the results obtained and the outline of the future iterations of 
assessments are given in this section. 

1 General description of the e^p elastic scattering 

The differential cross section of the elastic e^p scattering in the laboratory reference frame is 



parameterized by the well known Rosenbluth formula |11] 



da ( da 
dn ~ [dn 

where 



i + T 



4m2' \dn]^^ \2EJ sin4(0/2) 
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E is the incoming electron energy, 9 is the scattering angle of the electron in the lab frame, and 
nip is the proton mass. 

The Rosenbluth formula is based on the one photon exchange approximation QED diagram 
with modified proton-photon-proton vertex for the extended proton 

with Dirac (Fi) and Pauli {F2) FFs introduced. These FFs are connected with so called "Sachs" 
FFs as follows: 

In the non-relativistic limit they describe an electric charge and magnetic moment distributions 
within the proton with the normalizations as follows: 

G^(0) = 1, Gl,{0) = fip, 

G^(0) = 0, Gl,{0) = ixn: 
where Hp^n are the (p, n) magnetic moments and 
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is the mean square radius of the proton charge distribution in units [GeV ^j. In terms of 
invariant variables s and t the formula for differential cross section takes the form 



where s = s 



da 

2 

- mz., 



t tmi 



1 + - + 

s 



S2 



1 + r 



t = -q'^ = -AE'^ sm.\e/2) -[l + j^ sin2(6'/2)]"^ Note, that t varies in the 



+ 



p2 



2 s{s + + tml 



(2) 



range —s^/s < t < 0. When s — > oo, then one has 



da 



47rQ^ 
1^ 



1 + r 



(3) 



for any fixed value of t. It can be shown that the righthand side of (3) is an upper envelope 
curve for all 

As is seen from Fig. 1, this is confirmed experimentally in the wide ranges of both variables 
s and t. This is a direct evidence that one photon approximation works well at the available 
level of accuracy of the world experimental data. 

To model FFs, we follow a conventional way with introduction of the normalized to unity 
at t = isospin FFs, namely, -Fi'2 • They are connected with the Ff'2 and G^^\j through the 
relations: 



2 ) 



F? = \iFl-F^), F- = ^I^Fi 



f: 



2 ! 



where are the anomalous magtnetic moments of proton and neutron. 



2 Data sample 

We use the total sample of data on differential cross-sections for the electron-proton elastic 



scattering extracted from the PPDS ReacData data base |]2^ (see also Reaction Data [28| ). 
The complete list of references for the used experimental data is presented in Table 5.1. 

The data are divided by the sets corresponding to different experiments or experimental 
methodics. For example, if the different methodics were used in the same work to obtain the 
data, we treat different data samples as independent because of possible difference in systematic 
normalization errors. The data of sets 34 and 35 are taken from the graphs. In all the cases, 
where it was possible, we separated the systematic normalization errors from the total ones in 
order to renormalize the different data sets in a such a way that makes the whole compilation 
self-consistent. In a few cases we failed to separate the normalization errors from the total 
errors. We classify such data sets as "incomplete data" and will use them only at the stage of 
the first crude adjustments of the models. The incomplete data sets are sets 2,3,18,19,23,24,41 
of Table 5.1. They have only the total errors assigned without detailed information about 
systematic errors. 
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We have also excluded 5 points of set 1 (BUMILLER 61, PR 124, 1623, see Table 2.1.) that 
have inconsistent information about the errors, and few points that turned out to be in sharp 
contradictions with the nearest experimental points and are out of all parametrizations by more 
than 5-6 standard deviations. 



Table 2.1. Excluded points 



Experiment 




E[GeV] 


e\deg\ 


1 


BUMILLER 61 


0.76 


0.700 


145 






0.45 


0.800 


60 






0.33 


0.850 


45 






0.52 


0.875 


60 






0.70 


0.900 


75 


5 


YOUNT 62 


0.202 


0.307 


130 


6 


+YOUNT 62 


0.202 


0.307 


130 


25 


MISTRETTA 69 


0.422 


4.529 




26 


BREIDENBACH 70 




7. 


6 


37 


BORKOWSKY 74 


0.76[/m-^] 
0.56[/m-2] 


0.1498 
0.1801 





The resulting data sample that we have included in the first crude analysis (see section 5) consists 
of 847 experimental points covering the q range from 0.005 to 31.280 GeV'^ and 1961-1994 time 
period (see Fig. 1). 



3 Models 

To construct a simple but flexible enough estimator for the nucleon electromagnetic FFs, we 
explore the majority of previous old and recent ideas invented in the papers devoted to FFs 
modelling and data analyses [Q|,||2| - |24]: the simplest analytic structure of Pauli and Dirac 
FFs with the only established singularities that have a clear physical interpretation and with 
an asymptotic behaviour at a large momentum transfer as predicted by perturbative QCD [25|. 
Isoscalar and isovector FFsQ are modelled as follows: 

= i^GyAfD(i,2)(*)x^i;2W, (i) 

where F^y]y^£)(^i 2)(0 are the generalized vector dominance model factors, and Al'2 are the terms 
which control the asymptotic behavior of FFs at (— t) — > 00. 

All the functions used are the real analytic functions with real magnitudes at t < 0, with cuts 
dictated by the first few particles thresholds at t > 0. -^GV'Af_D{i 2)(*) have resonant poles only 
on secondary sheets of the Riemann surface of the FF and constructed as a sum of individual 
contributions corresponding to resonant poles. 

is the term motivated by possible poles of pionium spectrum (correspon- 
ding to 1 electromagnetic bound states of the tt+tt" system). There is 



1 ptomum 



(t) 



no possibility to get a parametrization of the whole pionium stectrum contribution to the FFs. 
Therefore, we introduce only the nearest to the point t = pole contribution to the FFs 

, _ (x/1 - t/4ml) • (1 + yi - (2m^ - e^)y4ml) 

Jpiomum{>^) — p. — n 9- , n IT. vTTl 9 ' ^ > 

V^l - tj^ml + x/1 - (2rn^ - ^^fl^ral 

It is the real analytic function in the whole t-plane with a righthand cut [4m^,oo[ with a pole 
on the second Rieman sheet. It corresponds to the lowest 1 pionium "bound" state with the 

'^Mathematica and Fortran codes of all the models can be obtained by request from: polislichuk@mx.ihep.su. 
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binding energy — ejr, which is unstable under annihilation. We neglect the pionium annihilation 
width in form (2). By construction it has the following properties: 

fpionium(fi) — 1) fpioniumii) ^ 1 at \t\ > CX), /pionmm (i) |e^=0 = !• 



is the model of p meson pole contribution to -^^yjvfD{i 2) whole t-plane. 

It is constructed from the real analytic function ^'^{t), proposed by Wataghin[p^]. It 
has the cut [4m^, cxd[ and two poles in the unphysical Riemann sheets in the complex conjugated 
points 

a>^(t) = "^^^ (3) 



where p pole positions in two complex conjugated points are both on the unphysical sheet as 
required by real analyticity. The pole positions are determined by requirement to reproduce the 
Breit-Wigner shape in the vicinity of p resonance {t = rn?^) at the upper edge of the cut in the 
small width approximation. The nonzero width p meson contributions to the -Pgv'MD(i 2)(^) 
given by 



with the following properties 



777 r 

r(0) = l, fP{t)^^ at |t|^oo. 



(t) ^'■^ models of oj meson contributions to FFs in the whole t-plane with the 

^'^ righthand cuts. They are the real analytic finctions with cuts [Am^,oo[ and 

[9m^,oo[ with poles on the unphysical Riemann sheets in the complex conjugated points, and 
are constructed as follows^ : 



$^2(t) = cosG^; 



p(3) 



ml-t + m^r^/vr • In ^T%^i 



r(2) 

sm 9^,2 ^, (5) 

where 00 resonance pole positions are in the complex conjugated points on the unphysical sheets 
corresponding to two and three body branch cuts. The pole positions are determined by the 
requirement to reproduce the Breit-Wigner shape in the vicinity of resonance on the upper edge 
of the cut in the small width approximation. The nonzero width omega meson contributions to 
the FFs are as follows: 



with properties 



/r:2(0) = l, /i"2(i)~^ at |t|-oo. 



^Here we try to test whether differrent branch channels can be identified from the fit. So, our form for uj and 
contributions differs from that used in p^,]2^,p4|. 
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fUt) 



are the models of (j) meson contributions to FFs in the whole t-plane with the 
righthand cuts. They are the real analytic finctions with cuts [4m|^,oo[and 
[9m^, (X){ with poles on the unphysical Riemann sheets in the complex conjugated 



points. They are constructed in a similar way 

^Ut) = cosef,- 



(3) 



2-t + m<^Vvr.lnJ^ 



m 



+ sin @l 2 



(2) 



t + m^T^ ■ 



m^.—Ami^ 



The nonzero width 6 meson contributions to the FFs have the form 



(7) 



with properties 



(8) 



const 
t 



at \t\ 



oo. 



Other ingredients of our hybrid model are multiplicative factors which control the FFs 
behaviour at large \t\. We use the two forms for it. 



are motivated by the parametrization of Mack [1 



1 - t/t 



.p 

cut 



1-t/t 



cut 



4mp is the proton-antiproton threshold, t 



cut 



'In 



9mi 



it 



s,v 
cut 



A2 



1,2, 



(9) 



''cut 



4m^, and A is the 



where 

scale parameter which is expected to be of the order of the Aqcd- 

M^'2{t) are the real analytic functions with singularities only on the righthand cuts [icnt;Co[ 
and [4mp,oo[, normalized by M^'2(0) = 1. At large —t they decrease faster than predicted by 
QCD, but they give quite good fits at available t as it will be shown further. This model is in- 
cluded in the current analysis in order to test whether the QCD asymptotic is really prominent 
at the available highest q'^. 

are motivated by parametrization of Gari and Kriimpelman [^] with taking 
into account the QCD asymptotic logarithms by analogy with the parametri- 



zations used in [20|,[23| 



(1 - t/fcut) 



i-vr+^rc^-t/t 



cut 



In 



"cut 



A2 



/ln%i 



2+u 



1,2, 



(10) 



GK^'"{t) are the real analytic finctions with singularities only on the righthand cuts [t^^j,oo[ 
and [4mp,oo[, normalized by Gi^r^'2(0) = 1 at t = 0. At large —t, they reproduce the FFs 
behaviour at large —t as predicted by pQCD. The exponent u is given |2H| by 



4 

3^' 



13=11 



-Nf 



0.174 for Nf = 5. 
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Now we are ready to assemble parametrizations for isotopis form factors. 
Modified Maclt-67 parametrization (Mack-98) 

= {{l-C¥,2)-fpionium{t)+Cf,2-flM'<Ml,2^ (H) 

with 10 adjustable parameters: Cf^'^, 0i',25 ^1^2- 

Modified Gari-Kriimpelman parametrization (GK-98-2V) 

-^1,2 ~ ((1 ~ ^1,2) ' fpioniumit) + Ci^2 " fl,2i^)) ^ ^^l,2i^)j 

FI2 = ((l-Cf^2)-/p— m(i) + Cf,2-r(*)) (12) 
with 10 adjustable parameters: Cf'^, Qi^2^ Vi%- 
Modified Gari-Kriimpelman parametrization (GK-98-3V)) 

-^1,2 = (i^ ~ C12 — 0^2) ■ fpionium{t) + C^2 ' /l^2(*) + C'f 2 • /i^2(0) ^ G'-^l,2(*)) 

Fl,2 = [{l-C',,2)-fpionium{t) + CP^2-nt))xGKl2 (13) 

with 14 adjustable parameters: Cf'^'"^, Oi^2j ®i 2' '7i'2- 

4 Preliminary fit results 

To select "the best models" for the FFs, we start with the crude fits of each above parametrization 
to the overall data sample (847 data points, see Table 5.1.). We used a mean least squares 
method with weights defined by the "total error" - the quadratically combined all kinds of 
experimental errors, ignoring correlations. The results of our best fits of analytic FFs are given 
in Table 4.1. The best fit estimations for the parameters and two kinds of their error estimates 
are presented: 1) the Minuit errors, 2) the errors obtained by propagation of the combined errors 
of the experimental data points to the parameter errors using equations to find the position of 
minimum in the parameter space. Since we will use the error propagation technique to get 
separate estimates of the parameter errors of different kinds due to statistical, systematic, and 
normalization errors in the input data, we must be sure that we are in a local minimum and 
parameters estimates are sufficiently close to the minimum position. These two errors must be 
close to each other if the usual Minuit assumptions about the minimum are valid and if the 
distance to the true minimum is sufficiently small (see Appendix 2). 
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Table 4.1. All data, weights made from all errors combined in quadratures. 



Model name 


Free & derived 


T — 1 I'll 1 1 111* 

Estimated parameter values and their errors 


& Fit quality 


parameters 




Mack-98 


xf 


2.702± 0.024(0.024) 




x\ 


0.894 ± 0.032(0.031) 






0.000 (fixed) 






0.05232 ± 0.00047(0.00047) 


9/7 r 


Cf 


0.9681 ± 0.0069(0.0068) 


986.4/836 


CP 


0.9127 ± 0.0049(0.0048) 


1.179 




0.5701 ± 0.0071(0.0071) 




6'"[deg] 


-0.75 ± 0.17(0.17) 




91.541± 0.069(0.079) 




6'2[deg] 


91.41413± 0.00013(0.00014) 




A[GeV] 


0.1024 ± 0.0039(0.0039) 




P rr 1 

r^[fmj 


0.860 ± 0.003(0.003) 


GK-98-2V 


Vi 


0.06699 ± 0.00142(0.00141) 




V2 


0.0364 ± 0.0011(0.0011) 




Vi 


0.1191 ± 0.0027(0.0027) 




V2 


0.0471 ± 0.0029(0.0029) 


X /dof 


C{ 


2.380 ± 0.041(0.041) 


973.6/836 


CP 


1.0843 ± 0.0093(0.0093) 


1.165 




-1.952 ± 0.038(0.038) 




C2 
fi [aegj 


8.01 ± 0.18(0.18) 




yi.o( 1 1 ± u.vvziiv.iiiiZL) 




^^[deg] 






Al vjc V J 


u.uioo HI u.uuou^u.uuou^ 




' E [^"^J 


n Sfi.'i + 002^0 002") 


GK-98-3V 


'/I 


06806 ± OOOQQfO 000991 




12 


0.03002 ± 0.00046(0.00047) 




1 L 


0.10627 ± 0.00048(0.00048) 






0.04876 ±0.00057(0.00057) 


X^/dof 




2.312 ± 0.026(0.023) 


972.4/833 


C2 


1.2383 ± 0.0070(0.0059) 


1.167 




-2.583 ± 0.023(0.021) 






8.99 ± 0.13(0.07) 




cf 


1.290 ± 0.022(0.022) 




^r[deg] 


1.91 ± 0.18(0.19) 




91.3779 ± 0.0029(0.0029) 




^t'[deg] 


91.40402 ± 0.00097(0.00097) 




<[deg] 


157.98 ± 0.99(0.99) 




^2*[deg] 


157.95 ± 8.15(8.11) 




A[GeV] 


0.0136 ± 0.0046(0.0046) 




r^[fm] 


0.861 ± 0.002(0.002) 


Polynomial 


ai 


-0.241 ± 0.005(0.005) 






0.023 ± 0.009(0.009) 


0.928 


r|[£m] 


0.849 ± 0.008(0.008) 



9 



In addition to the parametrizations under consideration we will also use the traditional 
second order polynomial 

1 + ai{- — 77) + a2 



4m2 4m2 

in the range from to 0.06 GeV^. Here 4m^ = 0.078 GeV^ is the bound of convergence circle 
for the presentation of FF by power series; ai and 02 are the free parameters. This parametriza- 
tion is often used to estimate the proton electric radius. It is sensitive to the behaviour of FFs 
near a zero value of and large deviation of "polynomial" value from predictions of other mod- 
els means that smooth parametrizations do not reproduce a possible fine structure near q^ = 
and probably may give wrong results for proton radius. On the other hand, such a "fine struc- 
ture" might be caused by inconsistency in normalizations of different data sets in this region. 
Such an inconsistency of the existing data in normalizations cannot be removed without using 
smooth parametrizations. Or we will have to wait until a new single experiment will precisely 
and densely enough measure the range of small g^, sufficient to get reliable polynomial fits. 

It should be noted that we obtain the best quality fits of the whole data sample of the 
da/dt by each of our three analytic models ever obtained by other analyses with other models. 
We quote in the tables the parameters errors estimated by two independent methods as it was 
described above. 

But all these fits are still not good enough to draw conclusions about the reliability of the 
physical parameters estimates (in spite of that we have obtained the estimates of the proton 
electric radius in a good agreement with the latest published estimates p3[). 

It should be stressed also that the model with three vector mesons has no advantages in 
describing the data in the space-like region of these naive fits. Pionium, p, and uj contributions 
are sufficient to get a relatively good data description. 

Let us try to get more reliable fits by excluding correlated normalization errors from the 
weights, because this may cause (as we will show further) a significant biases in parameter 
estimates. In the next section we use the correlated normalization uncertainties to make data 
sets relatively consistent with each other on the basis of the selected models. 



5 Data normalizing 

For further analyses we have excluded the so called "incomplete data sets": 2, 3, 18, 19, 23, 24, 
41 from the total sample of the data (see Table 5.1). They do not contain any information about 
systematics, and cannot be used in the procedure of renormalizations. To uncover a possible 
effect of this formal filtration we repeate the previous fit on the reduced sample of experimental 
data with the results presented in Table 5.2. As is seen we obtain quite good fit qualities for all 
the models and parameters shifted only slightly. As for the proton radii, they are not practically 
shifted. 

Having obtained a good enough fits by all the parameterizations, we remain on the firm 
ground to discuss the values of the proton radius. As is seen from Table 5.2, the model dependent 
radii are close to each other with overlapping error bands, but the radius obtained from a 
polynomial fit is smaller and its interval estimator, 

e [0.841, 0.857] [fm], 

does not overlap with a joined interval estimator of model radii, namely 

e [0.858, 0.867] [fm]. 
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Table 5.1. Statistics of data sets, renormalization factors, and references. 



Set 


Set 


After 


Normalization 


Points in 






No 


number of 


filtration 


factors 


90% CL 


REACDATA ShortCode 


Reference 




jjoints 






sample 






1 


53 


53 


0.985 ± 0.002 


45 


BUMILLER 61 


PR 124, 1623 


2 


9 





- 





OLSON 61 


PRL 6, 286 


3 


20 





- 





LITTAUER 61 


PRL 7, 141 


4 


2 


2 


0.999 ± 0.003 





LEHMANN 62 


PR 126, 1183 


5 


3 


3 


1.009 ± 0.004 


3 


YOUNT 62 


PR 128, 1842 


6 


3 


3 


1.000 ±0.000 


3 


+ YOUNT 62 


PR 128, 1842 


7 


8 


8 


1.017 ±0.004 


7 


DRICKEY 62 


PRL 9, 521 


8 


6 


6 


1.006 ± 0.001 


4 


DUDELZAK 63 


NC 28, 18 


9 


21 


21 


1.011 ± 0.011 


17 


BERKELMAN 63 


PR 130, 2061 


10 


7 


7 


0.988 ± 0.013 


5 


DUNNING 63 


PRL 10, 500 


11 


6 


6 


1.030 ±0.012 


5 


CHEN 65 


PR 141, 1267 


12 


12 


12 


1.052 ±0.009 


8 


+ CHEN 65 


PR 141, 1267 


13 


10 


10 


1.000 ± 0.000 


8 


FREREJACQUE 65 


PR 141, 1308 


14 


3 


3 


1.008 ± 0.006 


3 


GROSSETETE 65 


PR 141, 1435 


15 


93 


93 


1.000 ± 0.000 


88 


JANSSENS 66 


PR 142, 922 


16 


3 


3 


1.022 ± 0.005 





BENAKSAS 66 


PR 148, 1327 


17 


11 


11 


1.000 ± 0.000 


10 


BARTEL 66 


PRL 17, 608 


18 


10 





- 





ALBRECHT 66 


PRL 17, 1192 


19 


7 





— 





+ ALBRECHT 66 


PRL 17, 1192 


20 


24 


24 


0.986 ± 0.004 


20 


BEHREND 67 


NC A48, 140 


21 


4 


4 


0.912 ±0.013 


1 


± BEHREND 67 


NC A48, 140 


22 


1 


1 


1.009 ±0.015 


1 


BARTEL 67 


PL 25B, 236 


23 


4 





- 





BARTEL 67B 


PL 25B, 242 


24 


3 





- 





ALBRECHT 67 


PRL 18, 1014 


25 


4 


4 


0.997 ± 0.003 


3 


MISTRETTA 69 


PR 184, 1487 


26 


8 


8 


1.003 ± 0.009 


8 


BREIDENBACH 70 


MIT-2098-635 


27 


22 


22 


1.004 ± 0.003 


18 


LITT 70 


PL 31B, 40 


28 


15 


15 


1.002 ± 0.003 


11 


GOITEIN 70 


PR Dl, 2449 


29 


54 


54 


0.996 ± 0.002 


52 


BERGER 71 


PL 35B, 87 


30 


9 


9 


1.005 ±0.003 


9 


PRICE 71B 


PR D4, 45 


31 


4 


4 


0.985 ±0.012 


4 


GANICHOT 72 


NP A178, 545 


32 


5 


5 


1.003 ± 0.004 


5 


BARTEL 73 


NP B58, 429 


33 


16 


16 


0.997 ± 0.003 


15 


+ BARTEL 73 


NP B58, 429 


34 


37 


37 


0.992 ± 0.001 


35 


BOTTERILL 73C 


PL 46B, 125 


35 


30 


30 


0.988 ± 0.001 


22 


+ BOTTERILL 73C 


PL 46B, 125 


36 


11 


11 


0.997 ± 0.006 


10 


KIRK 73 


PR D8, 63 


37 


44 


44 


0.994 ± 0.001 


37 


BORKOWSKY 74 


NP A222, 269 


38 


11 


11 


1.013 ± 0.002 


11 


MURPHY 74 


PR C9, 2125 


39 


77 


77 


0.996 ± 0.001 


69 


BORKOWSKY 75 


NP B93, 461 


40 


7 


7 


0.993 ± 0.004 


7 


STEIN 75 


PR D12, 1884 


41 


3 





- 





MARTIN 76C 


PRL 38, 1320 


42 


8 


8 


1.022 ±0.008 


8 


MESTAYER 78 


SLAC-214 


43 


49 


49 


1.002 ±0.000 


40 


AKIMOV 78B 


YF 29, 922 


44 


5 


5 


1.002 ± 0.001 


5 


SIMON 80 


NP A333, 381 


45 


22 


22 


0.998 ± 0.002 


21 


WALKER 89 


PL 224B, 353 


46 


11 


11 


0.995 ± 0.005 


11 


BOSTED 89 


PR C42, 38 


47 


5 


5 


0.969 ± 0.011 


5 


ROCK 91 


PR D46, 24 


48 


13 


13 


0.994 ± 0.004 


12 


SILL 93 


PR D48, 29 


49 


22 


22 


0.998 ± 0.002 


22 


WALKER 94 


PR D49, 5671 


50 


32 


32 


1.000 ± 0.002 


28 


ANDIVAHIS 94 


PR D50, 5491 


All 


847 


791 




696 
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Table 5.2. Fits to sets of complete data only, weights made from all errors combined in quadratures. 



Model name 


Free & derived 


Estimated parameter values and their errors 


& Fit quality 


parameters 




Mack-98 


xl 


2.587 ± 0.091(0.079) 




x\ 


0.978 ± 0.072(0.061) 






0.0035 ± 0.0029(0.0027) 




X2 


0.0505 ± 0.0017(0.0015) 






0.9665 ± 0.0072(0.0071) 


X^/dof 


CP 


0.9063 ± 0.0051(0.0050) 


757.4/780 




0.5767 ± 0.0073(0.0072) 


0.971 




-0.97 ± 0.17(0.17) 




91.523 ± 0.045(0.044) 




^2[deg] 


91.414162 ± 0.000085(0.000084) 




A[GeV] 


0.1010 ± 0.0039(0.0039) 




r^[fm] 


0.861 ± 0.003(0.003) 


GK-98-2V 


rtl 


0.0624 ± 0.0015(0.0015) 




V2 


0.02762 ± 0.00090(0.00090) 




Vi 


0.1094 ± 0.0027(0.0027) 




Vl 


0.0427 ±0.0023(0.0023) 




Cf 


2.317 ± 0.038(0.039) 


750.4/780 




1.1894 ± 0.0090(0.0090) 


0.962 




-1.723 ± 0.035(0.035) 






9.01 ± 0.19(0.19) 




^i^[deg] 


91.3779 ± 0.0023(0.0023) 




^^[deg] 


91.40402 ± 0.00047(0.00047) 




A[GeV] 


0.0194 ± 0.0064(0.0064) 




r^[fm] 


0.865 ± 0.002(0.002) 


GK-98-3V 


Til 


0.0505 ± 0.0014(0.0014) 




Til 


0.02281 ± 0.00072(0.00073) 




Vt 


0.0989 ± 0.0020(0.0020) 




ri2 


0.0437 ±0.0017(0.0017) 


X^/dof 


C{ 


3.3392 ± 0.0092(0.0092) 


747.9/776 


CP 


1.3402 ± 0.0028(0.0028) 


0.963 


Cf 


-3.9995 ± 0.0092(0.0106) 




C^ 


9.000 ± 0.095(0.109) 




cf 


1.9176 ± 0.0061(0.0038) 






3.999 ± 0.056(0.034) 




91.37789 ± 0.00099 (0.00099) 




e^[deg] 


91.40399 ± 0.00051 (0.00051) 




<[dcg] 


157.99 ± 0.68(0.68) 




Oti^eg] 


158.0± 3.7 (3.7) 




A[GeV] 


0.0203 ± 0.0059(0.0060) 




r|[fm] 


0.862 ± 0.002(0.002) 


Polynomial 


ai 


-0.2407 ± 0.0047(0.0047) 


xVdof 


a2 


0.0229 ± 0.0088(0.0089) 


0.928 


r|[fm] 


0.849 ± 0.008(0.008) 
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Table 5.3. Fits to sets of complete data only, weights constracted without correlated normalization errors. 



Model name 


Free & derived 


Estimated parameter values and their errors 


& Fit quality 


parameters 




Mack-98 


xl 


2.03 ± 0.17(0.15) 




x\ 


1.89 ± 0.02(0.02) 






0.0755 ± 0.0055(0.0049) 






0.0225 ± 0.0015(0.0013) 






-0.2123 ± 0.0087(0.0086) 


1474.1/780 


CP 


0.6689 ± 0.0058(0.0057) 


1.889 




-1.699 ± 0.016(0.016) 






-9.37 ± 0.25(0.25) 




91.543 ± 0.017(0.017) 




^2[deg] 


91.414124 ± 0.000017(0.000017) 




A[GeV] 


0.00105 ± 0.00040(0.00036) 




r^[fm] 


0.909 ± 0.001(0.001) 


GK-98-2V 


rtl 


0.00589 ±0.00012(0.00007) 




V2 


0.004263 ± 0.000023(0.000041) 




Vi 


0.03401 ± 0.00064(0.00035) 






0.01211 ± 0.00043(0.00040) 






-10.47 ± 0.19(0.10) 


1447.7/780 




0.478 ± 0.021(0.019) 


1.856 




-9.3037 ± 0.0041(0.0060) 






1.567 ± 0.047(0.042) 




^i^[deg] 


91.504 ± 0.011(0.010) 




^^[deg] 


91.415123 ± 0.000002(0.000003) 




A[GeV] 


0.219000 ± 0.000012(0.000012) 




r^[fm] 


0.9115 ± 0.0018(0.0014) 


GK-98-3V 


Til 


0.005887 ± 0.000028(0.000028) 




Til 


0.004263 ± 0.000062(0.000063) 






0.03401 ± 0.00014(0.00014) 






0.012101 ± 0.00044(0.00045) 


X^/dof 


C{ 


10.47 ± 0.21(0.12) 


1447.7/776 


CP 


0.476 ± 0.034(0.032) 


1.866 


Cf 


-9.31 ± 0.27(0.16) 




c^ 


1.571 ± 0.074(0.072) 




cf 


0.0015 ± 0.0823(0.0525) 






-0.02 ± 0.27(0.26) 




91.50415 ± 0.00067 (0.00067) 




e^[deg] 


91.415123 ± 0.000001 (0.000001) 




<[dcg] 


157.9 ± 3.7(3.7) 




Oti^eg] 


158.0± 17.5 (17.9) 




A[GeV] 


0.2190 ± 0.0059(0.0060) 




r|[fm] 


0.911 ± 0.001(0.001) 


Polynomial 


ai 


-0.2542 ± 0.0030(0.0030) 


xVdof 


a2 


0.0372 ± 0.0050(0.0050) 


2.187 


r|[£m] 


0.873 ± 0.005(0.005) 
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Table 5.4. Fits with set of normalizations determined for each model. 



Model name 


Free and 




and 


derived 


Estimated parameter values and their errors 


fit quality 


parameters 




Mack-98 


xi 


1.94 ± 0.17(0.15) ± 0.06(norm) 




xl 


1.85 ± 0.020(0.019) ± 0.008(norm) 






0.0863 ± 0.0057(0.0052) ± 0.0023(norm) 






0.0256 ± 0.0014(0.0013) ± 0.0006(norm) 


X^dof 


C{ 


-0.2035 ± 0.0087(0.0086) ± 0.0026(norm) 


1195.6/780 




0.6660 ± 0.0058(0.0058) ± 0.0017(norm) 


1.533 




-1.744 ±0.015(0.015) ± 0.008(norm) 






-8.74 ± 0.25(0.25) ± 0.09(norm) 






91.519 ± 0.011(0.011) ± 0.009(norm) 




^2[deg] 


91.414116 ±0.000019(0.000019) ± 0.000009(norm) 




A[GeV] 


0.00229 ± 0.00067(0.00061) ± 0.00028(norm) 






0.8950 ± 0.0014(0.0013) ± 0.0009 (norm) 


Polynomial 


ai 


-0.2582 ± 0.0030(0.0030) ± 0.0560(norm) 


xVdof 


£12 


0.0505± 0.0050(0.0050) ± 0.0695(norm) 


1.603 




0.8794 ± 0.0051(0.0051) ± 0.0954 (norm) 


GK-98-2V 


Vi 


0.005646 ± 0.000028(0.000029) ± 0.000280(norm) 




V2 


0.004363 ± 0.000068(0.000072) ± 0.000185(norm) 




rit 


0.03277± 0.00013(0.00014) ± 0.00150(norm) 






0.01262 ± 0.00050(0.00054) ± 0.00155(norm) 


xVdof 


CI 


9.457 ± 0.074(0.080) ± 0.578(norm) 


1177.9/780 




0.4516 ± 0.0075(0.0080) ± 0.0284(norm) 


1.510 


ct 


-8.213 ± 0.079(0.084) ± 0.627(norm) 




^1 


1.617 ± 0.019(0.020) ± 0.074(norm) 




^r[deg] 


91.50423 ± 0.00079(0.00079) ± 0.00862(norm) 




«2[dcg] 


91.415123 ± 0.000001(0.000001) ± 0.000003(norm) 




A[GeV] 


0.219637 ± 0.000033(0.000033) ± 0.000500(norm) 




r|[fm] 


0.8986 ± 0.0012(0.0016) ± 0.0041 (norm) 


Polynomial 


ai 


-0.2598 ± 0.0030(0.0030) ± 0.0560(norm) 


xVdof 


02 


0.0531± 0.0050(0.0050) ± 0.0695(norm) 


1.605 


r^[fm] 


0.8822 ± 0.0051(0.0051) ± 0.0951 (norm) 


GK-98-3V 




0.005343 ± 0.000026(0.000026) ± 0.000021(norm) 






0.003972 ± 0.000044(0.000044) ± 0.000020(norm) 




vl 


0.03281 ± 0.00013(0.00013) ± O.OOOlO(norm) 




V2 


0.01221 ± 0.00030(0.00030) ± 0.00014(norm) 


X^/dof 


C{ 


9.285 ± 0.070(0.069) ± 0.045(norm) 


1176.1/776 


C2 


0.3355 ± 0.0070(0.0069) ± 0.0033(norm) 


1.516 




-9.573 ±0.074(0.073) ± 0.046(norm) 




C2 


1.778 ± 0.018(0.018) ±0.009(norm) 




cf 


1.3524 ± 0.0024(0.0024) ± 0.0019(norm) 




ct 


-1.255 ± 0.041(0.041) ± 0.015(noi-m) 






91.50414 ± 0.00066 (0.00066) ± 0.00057(norm) 




e2[dcg] 


91.415123 ± 0.000001 (0.000001) ± O.OOOOOl(norm) 




<[dcg] 


162.4313 ± 0.0011(0.0011) ± 0.0009(norm) 




^f[deg] 


162.370± 0.010(0.010) ± 0.007(norm) 




A[GeV] 


0.224153 ± 0.000056(0.000056) ±0.000023(norm) 




r^[fm] 


0.9029 ± 0.0009(0.0009) ± 0.0007 (norm) 


Polynomial 


ai 


-0.2607 ± 0.0030(0.0030) ± 0.0560(norm) 


xVdof 


Cl2 


0.0525± 0.0050(0.0050) ± 0.0695(norm) 


1.650 




0.8838 ± 0.0051(0.0051) ± 0.0949 (norm) 
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Table 5.5. Fits with common set of averaged normalizations. 





T?ypp ^7- Hprivprl 

J- J. v_. v_. oc 1 V c_i 


T'lstimfitprl ■n?ir?iTnptPv v^iliips ?inr] tVipir prvor^ 


Rr Pit mialitv 






Mack-98 




2.03 ± 0.14(0.10) ± 0.03(norm) 






1 866 ± 017(0 014') ± OOSfnorml 




Xn 


0.0819 ± 0.0046(0.0035) ±0.0016(norm) 




X2 


0243 ± OOllfO 0008) ± 0003fnorm) 


/do f 




-0 2069 ± 0087f0 0086) ± 0026fnorm~) 


1199 8/780 


^2 


6682 ± 0058f0 0058") ± 0017fnorm") 


1.538 


^1 


-1 737 ±0 OlSfO 015) ± OOSfnorm") 




^2 


O.OO ZC W.ZitJlW.Zi(J I ZC W-Wi/l llUilil 1 






91 520 ± OllfO Oil) ± 009('norm) 

(_/ j_ • ?_/ ^ 1 yj ■ XXI \j « X X f 1 « yj y.j *_y i xx^x xxx i 






91.414116 ±0.000019(0.000019) ± 0.000009(norm) 




A[GeV] 


00174 ± 00046(0 00035) ± 00017(norm) 

\J .\J\J ±. 1 X _1_ \J .\J\J\Jj:\J \\J .KJKJKJfJfJ J _1_ W . W W W X 1 \ XXV,»X XXX J 




' E I J 


0.8975 ± 0.0012(0.0010) ± 0.0008 (norm) 


GK-98-2V 


'/I 


005690 ± 000028(0 000029) ± 000023(norm) 

yj • \j \j 'J yj -y \j — i— yj •yjyjyjyj^tJ i — i— xj.yjyjyjyj^fjyi-L'^LLi-Li 




'12 


004326 ± 000064(0 000066) ± 000028(norm) 




'/I 


03299± 00013(0 00013) ± OOOll(norm) 

\J tXjxj ^ \j \j — 1 — \j X I \j • X \j t — 1 — \j • XXI XXV./X xxx 1 




'12 


01239 ± 00047(0 00049) ± 00021 (norm) 

yj a \j ajkj kJ — 1 — \J . yj yjyjAi i \yj i kJ j — i — yj tyjyjyj ^ -L. \ ±±\j±. xxx i 


Ido f 


^1 


9 53 ± 17(0 15) ± 08(norm) 

tj 'XJrj _i_ yj . ±. 1 \ w.X(_/i _i_ w • y xxv^x xxx i 


1179 5/780 


^2 


4575 ± 0078(0 0078) ± 0034(norm) 


1.512 


^1 


-8.29 ± 0.18(0.16) ± 0.09(norm) 






1.601 ± 0.024(0.023) ± O.Oll(norm) 






91.50425 ± 0.00075(0.00075) ± 0.00065(norm) 




6*^ [deg] 


91.415123 ± 0.000001(0.000001) ± 0.000001 (norm) 




A[GeV] 


0.21963 ± 0.00019(0.00017) ± 0.00008(norm) 




r|[fm] 


0.8999 ± 0.0014(0.0013) it 0.0006 (norm) 


GK-98-3V 




0.005314 ± 0.000025(0.000025) ± 0.000020(norm) 






0.003992 ± 0.000046(0.000046) ± 0.000021(norm) 






0.03267 ± 0.00012(0.00012) ± O.OOOlO(norm) 






0.01236 ± 0.00032(0.00032) ± 0.00015(norm) 


xVdof 




9.27 ± 0.13(0.25) ± 0.16(norm) 


1173.8/776 




0.334 ± 0.012(0.023) ± 0.013(norm) 


1.513 




-9.56 ±0.14(0.27) ± 0.17(norm) 




cf 


1.783 ± 0.029(0.051) ±0.031(norm) 
1.3524 ± 0.0024(0.0024) ± 0.0018(norm) 




^r[deg] 


-1.261 ± 0.042(0.042) ± 0.015(norm) 




91.50412 ± 0.00066 (0.00066) ± 0.00057(norm) 




^^[deg] 


91.415123 ± 0.000001 (0.000001) ± 0.000001 (norm) 




0f[deg] 


162.4313 ± 0.0011(0.0011) ± 0.0009(norm) 




0|[deg] 


162.371± 0.010(0.010) ± 0.007(norm) 




A[GeV] 


0.224163 ± 0.000056(0.000056) ±0.000023 (norm) 




r|[fm] 


0.9022 ± 0.0016(0.0030) ± 0.0021 (norm) 


Polynomial 


ai 


-0.2597 ± 0.0030(0.0030) ± 0.0560(norm) 


x'/dof 


02 


0.0517± 0.0050( 0.0050) ± 0.0695(norm) 


1.630 


r|[fm] 


0.8820 ± 0.0051(0.0051) ± 0.0951 (norm) 
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We treat this fact as a signal that our estimators are shghtly biased because all the fitted models 
are analytic with a true domain of analyticity and it has been expected that for all reliable fits 
the interval estimators of the model radii would belong to (or at least, well overlapping) the 
"polynomial radius" interval estimator. 

Now, to show that possible biases arc indeed presented in the above estimates, wc repeat all 
the fits on the complete data sets only, but with the weights constructed without admixture of 
the normalization errors as described below. The results are shown in Table 5.3. They exhibit 
significant shifts in the adjustable parameters and radii values, however, one should hardly trust 
these estimates because of a bad quality of fits. The bad fits quality is a signal that either data 
from different experiments arc inconsistent with each other in normalizations or the models 
under consideration have the limited areas of applicability. So, we will try to construct a self- 
consistent data sample guided by the models and normalization errors in the data. Then, we 
will try to find areas of application for each model. 

With each of the complete experimental data sets we associate the adjustable renormalization 
parameter 1/ A* to simulate a possible systematic rescaling of the data of different experiments 
with respect to each other. These normalization parameters are a subject of evaluation by 
minimizing a least squares functional + "penalty function" 



where i is the number of experiment. A* is the normalization for z'th experiment, y * is the 
vector contained the experimental differential cross section values of the i'th experiment, t * is 
the vector contained the calculated cross-section values, is the inverted error matrix for the 
data of the ?'th experiment, is the number of points at the z'th set. 

Aj was calculated analytically at each step of minimization. For details see Appendix 1. 

The second term in (2) is the "penalty function" increasing the value when the normal- 
izations come out from unity. 

Let us denote the total point-to-point normalization error quoted by experimentalists as 
where j is the point number in the ?'th data set. Analogously, let Si{j) be the total sys- 
tematic error from all other sources combined in quadratures. 

We apply the following procedure: 

• If systematic errors due to normalization were declared by authors, we use cri{j) = 

• If is not given, but is given, then o"j(j) = Sj(j), with motivation that in this 
case we have a reference number and we adopt that ^'j(j) < In such a way we will 
obtain the upper limit of the normalization systematic errors, since we cannot separate 
the non-normalization part of error. 

• As has been indicated earlier we do not analyze the data which contained only total errors 
without any information about their sources. Certainly, we lose part of statistics ( about 
7 %), however, we will obtain more reliable results with full control of parameter errors. 

So, at the end of the minimization we would have three sets of normalization factors, one 
for each model, and the same number of normalization correlation matrices. 

The results of the fits with renormalizations are presented in Table 5.4. We obtain the 
slightly better fits for each model but they all are still bad. We attribute this badness of fits 
either to the limited areas of the applicability of the models used or to the overestimation of the 
experimental errors in some experiments. 
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Let us stress that now we have three different estimates of the radius resulting from the 
polynomial fits, one for each set of the renormalization parameters. Now, we have the situation 
when the interval estimator of the radii derived from polynomial fits for each set of normalizations 
contains the corresponding interval estimator for the model radii (as it should be when there 
are no rapid oscillations near t = 0). 

As all the models give a rough equally "good" description, we cannot choose the "best" model 
and the corresponding set of normalizations among others. To proceed to the determination of 
the areas of applicability of the models selected we construct a common set of normalizations 
following the procedure described in the Appendix 1 ("weighted averaging of normalizations", 
see Table 5.1). We will use this set of normalizations in the subsequent iterations of the fits. 
Their error matrix^ will be used in the calculations of the parameter errors induced by our 
renormalization procedure. 

To check that we did not destroy drastically the previous results with three different sets of 
normalizations, we repeat all the fits with obtained common renormalization parameters. 

The results are presented in Table 5.5. As was expected we did get the "slightly worse fits" 
than those previously obtained with the individual sets of normalizations, but the whole picture 
of the results remains unchanged. All the three models give practically the same goodness of fit 
and closer to each other values of the proton electric radii. 



6 An attempt to get the confidential areas of applicability 

From previous sections we see that the "naive" (all errors quadratically combined) fits to the 
sample of complete data sets give pretty good chisquares, however, the fits to the same data 
sample with weights constructed from total errors without admixtures of correlated normal- 
ization errors result in unacceptably large chisquares. Using the renormalizations of the data 
admissable by the correlated normalization errors in the data we slightly improve the quality of 
the fits but they still remain unacceptable. This may be due to the following reasons or their 
combination: 

• each model has own and restricted areas of applicability which should be determined; 

• experimental data have hidden systematics not quoted by experimentalists. 

One possibility to get reliable parameter estimates is to use PDG scale factor prescription 
1 12 1 and to inflate parameter errors with scale- factor resulted from the fits (but not the errors 
due to the play with normalizations). But this way says nothing about the area of models 
applicability and may be relevant only if the areas of applicability are known. 

So, we have to determine the areas of applicability first or discover the subsample of data 
which has the overestimated experimental errors. The way that can help to identify the areas 
of applicability of each model is the following "model driven filtration procedure." [ 



The whole physical region of kinematic variables log q'^ versus cos 6^ is divided into bins 
and each parametrization is fitted to the data set using normalizations obtained above. 

The value of /^points is evaluated for each non-empty bin, where Npoints is the number 
of points per bin. The bin is declared as belonging to "the good description area" (GDA) 
for a particular model if this value is less than some threshold value. Otherwise, this bin 
is declared as belonging to "the bad description area" (BDA). The threshold normalized 



^ Final arrays of normalizations and their error matrices obtained for each model and their "weighted averages" 
can be obtained by contacting polishchuk@mx.ihep.su. 

*In these variables the distribution of points is sufficiently uniform. 
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value is chosen in such a way that the confidence level (CL) calculated for "good 
description" subset would exceed some threshold value (say 50% or 90%). The procedure 
starts from sufficiently large k = X^/-^6m ~ ^' ^^^^ iteration all the points belong 

the GDA. Then, we decrease n by small steps, calculate CL at each value of the k, and 
determine the corresponding temporal GDA subsets until the threshold value is reached. 

• The "good" and "bad" data subsets are formed from the bins classified in such a way. 
The "good" subset includes the points belonging to GDA of all models and, similarly, the 
"bad" subset contains the points belonging to BDA of all models. 

The results of our standard fits of all models on the "good at 90% CL" data sample are 
presented in Table 6.1. We see that rejection the "bad" data points not much affected on the 
values of the parameters for all models. The interval estimators are slightly biased but the joint 
interval estimate for the radius before filtration and after filtration are well overlapped. 

Before the model driven filtration, we have 

G [0.8912, 0.9021] [fm] three models, 
r| G [0.7818, 0.9822] [fm] polinomial. 



After the model driven filtration, we get 

r| G [0.8955,0.9059] [fm] 



G [0.7847,0.9885] [fm] 



three models, 
polinomial. 



From these interval estimators (constructed using a linear combination of the experimental 
and normalization errors), we see that the polynomial estimator is inconclusive because it gives 
too wide interval estimator. 

To show the stability of the radii under variation of the threshold CL value used in classifying 
data points as "good" and "bad" , the radii obtained at two values of the CL = 50 % and CL=90 % 
are tabulated below. 



Model 


50% CL, 

713 points 
in [fm] 


90% CL, 

696 points 
r|; in [fm] 


Mack-98 

GK-98-2V 
GK-98-3V 
Polynomial 


0.8963 ± 0.0017 

0.8976 ± 0.0011 
0.9008 ± 0.0011 
0.8842 ± 0.0058 


0.8937 ± 0.0015 
0.8971 ± 0.0016 
0.9005 ± 0.0010 
0.8866 ± 0.0062 



The quoted errors do not include the contributions due to normalization errors. 

Now, let us discuss further the areas of applicability. If the subset of bad data is not empty 
and is grouped in such a way that this subset can be interpreted as the area where all models 
failed then the good data subset belongs to the area of applicability with a preselected CL. 

If the bad data points are randomly distributed over the whole available kinematic region, 
then we have only one possibility: To declare that the models are applicable in the whole region 
and the bad data are the data with overestimated errors (occasionally outstanding data). In 
this case, in order to estimate the errors in the parameters, we should use the Birge scale- 
factor method. The situation with bad data bins is presented in Fig. 3. There are no prominent 
regularities in the population of bad data bins in the whole kinematic region and we conclude that 
our models are applicable in the whole kinematic region, except for a few bins with presumably 
overestimated experimental errors. The only way to get reliable errors is the usage of the scale 
factor prescription. As the best estimate of the proton electric radius in terms of average, 
dispersion and systematic errors, we take the estimate derived from the GK-98-2V fit, taking 
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into account the scale factor \/1.512 (see Table 5.5), and assign the additional systematic error 
due to model selection as one half of the maximal difference between central values of the radii. 
So, we have 

= 0.900 ± 0.002(exp) ± O.OOl(norm) ± 0m2{models) [fm]. 

The same procedure of constructing a final radius estimate based on the filtered data (see 
Table 6.1) gives 

= 0.897 ± 0.002(exp) ± 0.001 (norm) ± 0m3{models) [fm]. 

Wc sec that the estimate obtained with the model driven filtration gives the result consistent 
with that obtained by using the Birgc factor. So, our best estimate of the proton radius is the 
estimate based on the model driven filtration. 
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Model name 


Free & derived 


Estimated parameter values and their errors 


& Fit quality 


parameters 




Mack-98 


xl 


1.59 ± 0.17(0.13)± 0.039(norm) 






1.806 ± 0.020(0.017)± 0.0056(norm) 




X2 


0.0953 ± 0.0057(0.0045)± 0.0018(norm) 




X2 


0.0284 ± 0.0014(0.0010) ± 0.0004(norm) 


xVdof 


Cf 


-0.213 ± 0.011(0.011) ± 0.003(norm) 


510.3/685 


CP 


0.6722 ± 0.0078(0.0077) ± 0.0019(norm) 


0.745 


ct 


-1.725 ± 0.020(0.020) ± 0.007(norm) 






-9.07 ± 0.30(0.30) ± 0.09 (norm) 






91.521±0.014(0.014) ±0.009(norm) 




^^[deg] 


91.414116±0.000021(0.000021) ±0.000008(norm) 




A[GeV] 
r^[fm] 


0.0035 ± 0.0009(0.0007) ± 0.0003(norm) 




0.8937 ± 0.0015(0.0013) ± 0.0008 (norm) 


GK-98-2V 


nl 


0.005652 ±0.000032(0.000032) ± 0.000021 (norm) 






0.004407 ± 0.000077(0.000077) ± 0.000027(norm) 




rii 


0.03275 ± 0.00015(0.00015) ± O.OOOll(norm) 






0.01253 ± 0.00055(0.00056) ± 0.00019(norm) 


xVdof 




9.25 ± 0.19(0.15) ± 0.08(norm) 


493.9/685 




0.457 ± 0.010(0.010) ± 0.003(norm) 


0.721 


ct 


-7.98 ± 0.20(0.17) ± 0.08(norm) 






1.594 ± 0.030(0.030) ± O.Oll(norm) 






91.50427 ± 0.00094(0.00095) ± 0.00068(norm) 




^^[deg] 


91.415123 ± 0.000001(0.000001) ± 0.000001 (norm) 




A[GeV] 
r%[im\ 


0.21991 ± 0.00023(0.00019) ± 0.00008(norm) 




0.8971 ± 0.0016(0.0014) ± 0.0008 (norm) 


GK-98-3V 


ril 


0.005844 ± 0.000032(0.000030) ± 0.000021 (norm) 




ril 


0.004495 ± 0.000073(0.000072) ± 0.000024(norm) 




Vl 


0.03358 ± 0.00015(0.00014) ± O.OOOlO(norm) 




Vl 


0.01253 ± 0.00050(0.00050) ± 0.00016(norm) 


xVdof 


Cf 


9.877 ± 0.077(0.072) ± 0.036(norm) 


494.8/681 




0.4161 ± 0.0089(0.0084) ± 0.0028(norm) 


0.726 


Cf 


-8.798 ± 0.082(0.076) ± 0.036(norm) 






l.Doo ± [).[)ZZ[{).[jZi) ± (J.L)L)y(normj 




Cf 


0.1937 ± 0.0053(0.0053) ± 0.0034(norm) 




ct 

0i[deg] 


-0.137 ± 0.072(0.072) ± 0.020(norm) 




91.50409 ± 0.00085(0.00085) ± 0.00062(norm) 




O2M 


91.415123 ± 0.000001 (0.000001) ± 0.000001 (norm) 




<[deg] 


157.8 ± 2.6(2.6) ± 1.9(norm) 




etideg] 


159.9± 13.0(13.0) ± 8.0(norm) 




A[GeV] 
r^[fm] 


0.219171 ± 0.000084(0.000085) ± 0.000038(norm) 




0.9005 ± 0.0010(0.0009) ± 0.0006 (norm) 


Polynomial 


Oi 


-0.2624 ± 0.0036(0.0036) ± 0.0567(norm) 


X^/dof 


0-2 


0.0575 ± 0.0059(0.0059) ± 0.0697(norm) 


0.817 




0.8866 ± 0.0062(0.0062) ± 0.0957(norm) 



Table 6.1. Results of the last iteration of fits on filtered data sample. 



7 Summary 

We have made the comparison of the majority of pubhshed parametrizations of the proton 
electromagnetic FFs using the largest experimental data sample on da/dt ever used. All the 
models give Tinacceptable "naive" fit quality to the whole data sample and we have to select and 
modify three parametrizations for further detailed analyses. These three hybrid parametrizations 
give the best and close to each other values of the x2/cZ.o./. parameter which characterizes the 
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fit quality (for the "naive" fits it is less than 1.2 for each model, see Table 4.1). 

All the reconstructed models have correct analytical structures with the minimal set of 
known singularities. They are all easily generalized for inclusion of the additional Wataghin 
terms |15|,|22| into the GVMD multiplier of the FF models for additional vector mesons as is 
seen from inspection formulas for the GK-98-2V and GK-98-3V parametrizations. It should be 
noted that the addition of i;^-meson did not resulted in improving the fits quality to the data 
in the space-like region (remember, that the addition of (/>-meson gives four extra adjustable 
parameters), but shifted slightly the proton radius towards higher values. 

The common feature of the tree models used is the presence of pionium contribution. In our 
preparatory attempts to produce cross assessments with models without these contributions, we 
had to reject by all the filtrations more experimental data than in selected and modified models 
with pionium. In the case with pionium all the filtrations rejected less than 18 %. It is strong 
evidence in favour of the pionium contribution to the GVDM multiplier of the FFs. The exact 
form of pionium contributions to the FFs is the task for the future refinements of our procedure 
and models. 

It should be stressed that in the recent experiment conducted by L.L.Nemenov with co- 
workers at IHEP (Protvino) the pionium was detected experimentally p9|]. 
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As it is seen from Tables 5.5 and 6.1 all three models give statistically equivalent descriptions 
of the unfiltered and filtered data. It can be also seen in Figs. 4-5, where the modulae of the 
electric and magnetic formfactors for each model are presented for comparison in a large interval 
of negative and positive q^. 

The model GK-98-2V is preferable at this stage of data and models cross assessment because 
it has a proper behaviour at high and gives the best description of the assessed data sample. 
The main results are summarized as follows: 

• We have developed the data and models cross assessments ( CRAS) technology in the case 
of the data on elastic e^p scattering and models of the nucleon form factors. We believe 
that the CRAS technology will be applicable in other well and will be helpful in 
the following tasks: 

— to make data from different experiments mutually consistent on the basis of available 
relevant models; 

— to determine the areas of applicability of the competitive models; 

— to produce technical and derived physical parameters estimations with metrological 
quality. 

• Using the CRAS technology we have constructed, from the well known inputs, the currently 
best parametrization of the nucleon FFs in the space-like region which is extendable in 
the sense of GVDM, can be also used to describe nucleon FFs in the time-like region, and 
is capable to eliminate the "false" contradiction between the charge radius of the proton 
from the hydrogen Lamb shift measurements (see recent review of S.G. Karshenboim [^) 
and from the high-energy electron-proton elastic scattering data. 

From the fit results of all the models to the "90% CL" sample, we have derived our best 
estimate of the proton charge radius 

r|; = 0.897 ± 0.002(exp) ± O.OOl(norm) ± 0.003(modeZs) [fm] 

given by the modified Gary-Kriimpelman model on the filtered data sample. This estimate is 
in good agreement with the estimate derived from the fits to unfiltered data but with the usage 
of the Birge method to justify the errors assigned to the adjustable parameters 

= 0.900 ± 0.002(ea;p) ± 0.001 (norm) ± 0.002(mode/s) [fm]. 

The traditional method of extracting proton radius from quadratic polynomial fit gives the 
inconclusive estimates on the model driven renormalized data (see Table 5.5 and Table 6.1), 
because they have too large systematic errors due to more correct treatment of the systematic 
error in the radius induced by the correlated experimental errors of normalization. 

It should be noted that the above results on the parameters of the three models analysed are 
still preliminary. They are a subject for the re-estimation with the data on n e~, n elastic 
scattering and e^e~ nucleon anti-nucleon added to the data on p elastic scattering. It is 
likely that all models will require the inclusion of the protonium atom singularities. We plan to 
perform further CRAS-iterations for the nucleon form factors with an extended data sample in 
the near future. 

A straightforward outcome of our analyses for the current and planned experiments are as 
follows: 

• The HERA data on e^p — > e^p would be extremely useful not only for the extraction of 
the FFs but also to check the relevancy of one photon approximation in the description of 
the da/dt data at moderate —t. Unfortunately, these measurements are not discussed or 
planned in the HERA physics program. 
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• The CEBAF data near t = will be very helpful to test the GVDM with pionium contri- 
bution (note that "bad" data points are concentrated near t = 0, see Fig.3.) 

To simplify the reproduction of our results, we prepare a special "WEB e^p — > e^p data site" 
with computer readable files of all data, models, and current results for the arrays of parameter 
values and their correlators on the basis of PPDS data. It can be reached by 

\protect\vrule widthOpt\protect\href {http: //sirius . ihep. su\string~polishch/ep 

or 

\protect\vrule widthOpt\protect\href {http: //dbserv. ihep. su\string~bvp/ep2ep_t 
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A Normalizations and their correlations 



We shall use the following notation: 

y ^ — array of average values of an observable y{x) measured at ra* x-points in i-th experiment. 

— X error matrix of z-th experiment. 

M^* = (T*)""*" — weight matrix for the data of i-th experiment. 

^i(j ) — estimation of normalization uncertainty in j-th point of the i-th experiment. 

A' auxiliary normalization factor to be assigned to i-th experiment to tune mutual normal- 
izations of different independent experiments. 

We shall use the following expression for x2 

i=i i=i j=i 

which must be minimized with respect to model parameters p and with respect to auxiliary 
normalization factors. Necessary requirements for extremums lead to equations 

= -2 E ^'"'^^ • ^'^y ' - ^'^'t'^i^) = 0, (2) 



- = -2f -(i) • W\y ' - A-'^f -(i)) + 2 V —-^ = 0. (3) 

1^1 nKj) 



Prom the second set of equations one can express auxiliary normalizations for m-th model as 
follows: 

Inserting these expressions for normalization factors into above expression for x2™, we will get 
a new formula for x2™, which must be minimized with respect to model parameters only. After 
solving the minimization task and selecting proper local minimum point (p™^), we can obtain 
the corresponding error matrix 

1 /^v9"* 

"dim "13m 

In a general case the normalization parameters for different models are correlated, and some of 
the normalizations may be strongly correlated for one model. Let us first prepare formulae for 
correlators between normalizations for different models. 

am Hi Pm' ^^rn' M,' 
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2 



^ aA-'^ aA-'-^- . , .2 , e V ^2 .o^ 

am/ij ^"^"^ C'yMj '^i'Mi "^^W 



where 



•q.v i^in ^P-i ^t^i 

and partial derivatives are estimated at minima points for m and m' models. Now we are able 
to construct estimators for a unique set of lambdas and their error matrix 

(9) 

and 

JV-. = rV— — r^fV— — V (lo) 

V ^ / J nimm ' V / j J\jmm ' / j ivrrnrn ivrm'm' ' ^ ' 

m m 3j mm' n 



B Error propagation 

B.l Errors of model pcirameters 

Having numerical values for normalization correction factors, we can find a final parameter 
estimation. We will use the standard expression for modified x2"^ 

TV 

X2™ = ^(y ^ - A"^''f'"(i)) • W\y'- A"^'*f '"(i)), 

i=l 

with axillary correction factors fixed at their best values obtained above. After minimization we 
will obtain new estimations for model parameters and corresponding parameter errors 

To estimate the contribution to parameters covariance from renormalization, we will propagate 

covariances NJ^ as follows: 

''J 

-f-rm,norm _ / \norm _ \ ^ '^Poim Arrra^^/3m 

Following the same way, we can propagate systematic (or statistical) point-to-point errors to be 
the errors of parameters: 

jTm,sys _ /A m A m \sys _ \ " ^^am ATt-'^^J Prn 
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here //j and Vj are the yu'th point in the i'th experiment and z^'th point in the j'th experiment, 
respectively. If usual MINUIT assumptions about the minimum are valid, this matrix must 
coincide with the MINUIT error matrix. Really, we have at the minimum position 

2 dpaopis 

in the absence of point-to-point correlation 



it corresponds to the standard MINUIT error definition. 

To calculate partial derivatives of parameters over lambdas, we will use the usual trick^ pl| . 
The set of equations 

g 2™ 
-TT — = 

might be considered as the set of functional equations defining dependencies of parameters upon 
lambdas. Taking the total derivative from both parts of above equations over A"*'*, we will 
obtain a linear system 



from which 



Finally, we will have 



jjm,norm _ _ \ ^ \ ^ -Qm A^"^ 11"^ 



1/ 



with 

92 x2'" 



dX'^'^dpr^ dp^ 



W\2X'^^H'^{%)-y'). 



Thus, for each model t*" we will have full estimations of parameters and their covariances 

/ m\ T-rm -r-rm,norm 



""We thank our colleague S.I. Alekhin who point out this method 
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B.2 Errors of proton cheirge radius 

Errors of the proton charge radius are calculated by propagation of the parameter errors of each 
type separately. The "systematic" error due to the uncertainties in normalizations 



TP -^^nrn rtrtrrr} • JP 

n 



m,{norm) ~ 2-^ Q rn oimPm Qjfn. 



The error due to the "total experimental" parameters errors obtained from MINUIT system 

a2 _ \ ^ _E_TTm '-" E 
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Figure 1: Data atlas. 
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Proton Charge Radius (fm) 
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Figure 2: The proton charge radius as it is usually known. (The original picture was in the 



article of S.S. Karshenboim, see hep-ph / 9712347 ) 
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Scattering angle distributiorTI 



Figure 3a: Populati 
as a function of 



atito of tlie "bad" data bins 
electron scattering angle. 
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Figure 3b: Population of the "bad" data bins 
as a function of momentum transferred. 
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Figure 3: Population of the bad data bins in the whole kinematic region. 
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